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ABSTRACT 

We present CARMA observations of the thermal dust emission from the circumstellar disks 
around the young stars RY Tau and DG Tau at wavelengths of 1.3 mm and 2.8 mm. The angular 
resolution of the maps is as high as 0.15", or 20 AU at the distance of the Taurus cloud, which 
is a factor of 2 higher than has been achieved to date at these wavelengths. The unprecedented 
detail of the resulting disk images enables us to address three important questions related to the 
formation of planets. (1) What is the radial distribution of the circumstellar dust? (2) Does the 
dust emission show any indication of gaps that might signify the presence of (proto-)planets? (3) 
Do the dust properties depend on the orbital radius? 

We find that modeling the disk surface density in terms of either a classical power law or 
the similarity solution for viscous disk evolution, reproduces the observations well. Both models 
constrain the surface density between 15 and 50 AU to within 30% for a given dust opacity. 
Outside this range, the densities inferred from the two models differ by almost an order of 
magnitude. The 1.3 mm image from RY Tau shows two peaks separated by 0.2" with a decline in 
the dust emission toward the stellar position, which is significant at about 2-Aa. For both RY Tau 
and DG Tau, the dust emission at radii larger than 15 AU displays no significant deviation from 
an unperturbed viscous disk model. In particular, no radial gaps in the dust distribution are 
detected. Under reasonable assumptions, we exclude the presence of planets more massive than 
5 Jupiter masses orbiting either star at distances between about 10 and 60 AU, unless such a 
planet is so young that there has been insufficient time to open a gap in the disk surface density. 
The radial variation of the dust opacity slope, /?, was investigated by comparing the 1.3 mm 
and 2.8 mm observations. We find mean values of (3 of 0.5 and 0.7 for DG Tau and RY Tau 
respectively. Variations in /3 are smaller than A/3 = 0.7 between 20 and 70 AU. These results 
confirm that the circumstellar dust throughout these disks differs significantly from dust in the 
interstellar medium. 



Subject headings: 

1. Introduction 

Resolved images of circumstellar disks around 
young stars provide the most direct tool for inves- 
tigating the formation of planets. At millimeter 
wavelengths, the thermal dust emission is gener- 
ally optically thin and meas ures the radial distri- 



butio n of circumstellar dust ([Beckwith fc Sargent 
Il990l) . However, since circumstellar disks in 



nearby star forming regions typically have radii 
between 100 and 500 AU, sub-arcsecond angu- 
lar resolution is required to spatially resolve the 
dust emission, even in nearby star-forming clouds. 
Millimeter-wave interferometers are essential for 
such studies. 

Since sub-arcsecond observations at millimeter 
wavelengths require both high sensitivity and high 
dynamical range, only a small number of bright 
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disks have been observed at resolutions of 0.4" 
1 



to date (Brown et al. 2008 



Guillote au et aL 
19991: llsella et al. Il2007t iPietu et al. 112005. 200 



20001 ). The Combined Array for Research in 
Millimeter-wave Astronomy (CARMA) and the 
new extended configuration of the Sub-Millimeter 
Array are rapidly enabling more extensive high 
resolution surveys of circumstellar disks, particu- 
larly i n the Taurus and Ophiuchus star forming re- 



gions (lAndrews et al. 2009t iHughes et al. 
llsella et al. Il2009l hereafter Paper I). 



2009t 



The highest angular resolution achieved so 
far by millimeter-wave interferometers is 0.3"- 
0.4", corresponding to spatial scales of 40-50 AU 
at the distance of Taurus and Ophiuchus. In 
most cases, the dust density appears to in- 
crease smoothly inward down to the orbital 
radius resolved by the observations, typically 
^^25 AU. However, central cavities in the dust 
distribution are revealed in a number of disks 
(|Andrews et al. I l2009t 



Hughes et aU l2009l ). It 



remains a matter of debate whether these cavi- 
ties are caused by dynamical interactions, inside 
out disk dispersal mechanisms, d ust opacity varia- 
tions , or visco us evolution (e.g., Alexander et aL 
20061: ICalvet ] |2005: Chia ng &: Murrav-Clav ..200 ' 



Dullemond fc Dominik ]l2005|. Paper I) 



a- 



Nevertheless, these observations still lack the 
angular resolution required to resolve the inner- 
most part of the disk where the density of the 
circumstellar material is highest and the forma- 
tion of planets is more probable. Here we de- 
scribe CARMA observations of the thermal dust 
emission towards the young stars DC Tau and 
RY Tau at an angular resolution of 0.15" at 
1.3 mm and 0.3" at 2.8 mm. At the distance 
of Taurus (140 pc), 0.15" corresponds to spatial 
scales of 20 AU, such that emission on orbital 
scales comparable to Saturn can be resolved. This 
is more than a factor of two improvement over pre- 
vious observations of circumstellar disks at these 
wavelengths. 

DC Tau and RY Tau are classical T Tauri 
stars of spectral t ype M O and Kl respectively 



(iMuzerolle et al. 1998t Kenvon fc Hartmann 



19951 ) . Stellar ages inferred from stellar evolution- 



ary models are less than 1 Myr (see Paper I for 
more details and references). The relative youth 
of both systems is confirmed by the presence of 



large amounts of gas and dust extending to 0.1 pc 
and by associated ste llar jets and outflows (see 



2007 : ISimon et al. koOollTesti et al. II 20031: fWiluCT et al. I \20m 



e.g.lMcGroartv fc Ra v 2004; St-Onge fc Bastien 
2008r). From near- infrared to millimeter wave- 



lengths, both objects exhibit strong emission in ex- 
cess of that from the stellar photospheres. This is 
attributed to rotating disks with radii of few hun- 
dred AU that first absorb and then re-emit radia- 



tion from the central s tars. (jKoerner fc Sargent 



19951: iTesti et al. Il2002[ ). Our earlier CARMA ob- 



servations of 1.3 mm thermal dust emission from 
these disks, at a resolution of 0.7", suggested disk 
masses between 5 and 150% of the stellar mass 
for both sources (see Paper I). These high disk 
masses and the youth of RY Tau and DC Tau 
make these prime targets to investigate the earli- 
est stages of planet formation. Our new observa- 
tions of RY Tau and DC Tau have a factor of 5 
better angular resolution and a factor of 3 better 
sensitivity than the previous data. 

This paper investigates three main questions re- 
lated to the formation of planets in young circum- 
stellar disks. (1) What is the surface density dis- 
tribution in the observed disks down to an orbital 
radius of 10 AU? (2) Are there any signatures of 
planet formation contained in the dust distribu- 
tion? Finally, (3), do the dust properties vary with 
orbital radius? A qualitative answer to the first 
two questions is proposed in Section [3] where we 
present the observations and discuss the morphol- 
ogy of the dust emission. A quantitative analysis 
is described in Section |4l where we compare the 
observations with theoretical models of disk emis- 
sion. Implications of these results for disk struc- 
ture, for the possible presence of planets, and for 
the radial variation of the dust opacity are consid- 
ered in Section [5] The conclusions are presented 
in Section [6] 

2. Observations and data reduction 

We observed thermal dust emission from the 
RY Tau and DC Tau circumstellar disks using 
CARMA in the A, B, and C configurations. The 
date of observation, array configurations used, 
baseline range, sizes and orientations of the syn- 
thesized beams, integrated fiuxes, seeing and 
noise levels are summarized in Table [TJ The C- 
configuration observations were presented in Pa- 
per I. 
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Table 1 

Summary of CARMA continuum observations 



Object Date Array Baseline Beam Flux Seeing Noise 

(UTC) Configuration range (m) FWHM("), PA(°) (mJy) (") (mJy/beam) 



Observations at 1.3 nun 



DG Tau 




Ah 


hBH 


-C 




0.17x0.15, 103 


367±14 0.03 


0.96 




2007 Oct 08 


c 






21-279 


1.04x0.81, 112 








2007 Dec 14 


B 






81-937 


0.33x0.29, 116 








2009 Jan 31 


A 






130-1884 


0.15x0.13, 102 






RY Tau 




Ah 


hBH 


-C 




0.17x0.14, 81 


227±7 


0.90 




2007 Oct 22 


c 






16-280 


1.24x0.78, 102 








2008 Dec 30 


B 






82-935 


0.40x0.31, 109 








2009 Jan 19 


A 






139-1884 


0.15x0.13, 82 






Observation at 2.7 mm 


DG Tau 




Ah 


KB 






0.45x0.38, 131 


58±6 0.07 


0.45 




2008 Jan 14 


B 






80-798 


1.12x0.60, 125 








2009 Feb 05 


A 






136-1678 


0.34x0.32, 163 






RY Tau 




Ah 


hB 






0.36x0.30, 82 


36±3 0.07 


0.28 




2008 Feb 01 


B 






82-945 


0.77x0.59, 95 








2009 Feb 10 


A 






123-1884 


0.35x0.29, 82 







The observations were obtained at LO frequen- 
cies of 228.1 GHz (A = 1.3 mm) and 106.2 GHz 
(A = 2.8 mm). The CARMA correlator at the 
time of the observations contained three bands, 
each of which was configured to 468 MHz band- 
width to provide maximum continuum sensitivity. 
The band pass shape was calibrated by observ- 
ing 3C273; flux calibration was set by observing 
Uranus and 3C84. The radio galaxy 3C111 was ob- 
served every 9 minutes to correct for atmospheric 
and instrumental effects. Variations of the atmo- 
spheric conditions on time scales shorter than 9 
minutes are not corrected, in effect resulting in 
seeing. We quantifled the atmospheric seeing by 
measuring the size of the phase calibrator image; 
if the seeing is negligible the phase calibrator ap- 
pears as a point source. Otherwise, the seeing pro- 
duces a Gaussian smoothing that can be quantified 
through the full width half maximum (FWHM) of 
the resulting image. We find that at 1.3 mm the 
effect of seeing is negligible for RY Tau but pro- 
duces a FWHM of 0.03" for DG Tau. Atmospheric 
conditions were slightly worse during the 2.8 mm 



observations, resulting in seeing of 0.07" for both 
objects. These seeing estimates do not account for 
variations in the atmospheric conditions on angu- 
lar scales of 10°, corresponding to the separation 
between the source and the calibrator. Values for 
the atmospheric seeing are summarized in Tabled] 
and are adopted in the model fitting described in 
Section m 

The raw data were reduced using the MIRIAD 
software package. The maps of the continuum 
emission shown in Figure [1] were derived using 
GILD AS software. Corresponding complex visi- 
bilities are shown in Figure [21 At 1.3 mm, natural 
weighting of the A, B and C configuration obser- 
vations produced a FWHM synthesized beam size 
of ~ 0.15". The noise levels are 0.96 mJy/beam 
and 0.90 mJy/beam respectively for DG Tau and 
RY Tau. Dust emission at 2.8 mm was observed in 
the A and B configurations at angular resolution 
of ~ 0.35" and noise levels of 0.45 mJy/beam and 
0.28 mJy/beam for DG Tau and RY Tau respec- 
tively. 
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6R.A. (arcsec) 6R.A. (arcsec) 

Fig. 1. — : Maps of the dust thermal emission observed at a wavelength of 1.3 mm (upper panels) and 2.8 
mm (lower panels) towards RY Tau (left panels) and DG Tau (right panels). The color scale shows the 
surface brightness starting from the 3cr level, with contours plotted every 4cr. The la noise level and the size 
of the synthesized beam are given in Table [TJ The inset in the upper left panel shows the central 0.4" x 0.4" 
region of the RY Tau disk where contours start at 28ct with increments of la. The surface brightness is 
characterized by two peaks separated by ^-^0.2". 



3. Morphology of the dust emission 

In Figure [U the dust emission in both disks 
is clearly resolved and characterized by a smooth 
and centrally symmetric radial profile. DG Tau 
intensity contours are almost circular suggesting 
a disk inclination smaller than 30°. For RY Tau, 
the intensity contours are elongated in the North- 
East direction suggesting a disk position angle of 



about 24° measured East from North and a disk 
inclination of at least 65°. For both sources the 
disk orientations agree with those found in Paper 
I. 

3.1. RY Tau disk morphology 

The 1.3 mm dust continuum emission from the 
RY Tau disk shows two spatially resolved peaks 
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Fig. 2. — : Correlated flux measured at 1.3 mm and 2.8 mm toward RY Tau (left panels) and DG Tau (right 
panels) as ftmction of the baseline length Buv expressed in kA. 



separated by about 0.2" (28 AU), and oriented 
along the apparent major axis of the disk. De- 
tails of the central 0.4" x 0.4" region are displayed 
in the inset in the upper left panel of Figure [1] 
and the radial profile of the surface density along 
the disk major axis is shown in Figure [3] The 
intensity at both peaks is 29 mJy/beam, which 
is 2 mJy/beam (i.e. 2.2(t) higher than the inten- 
sity at the center of the disk. We also estimated 
the expected central surface brightness by fitting 
a gaussian to the surface brightness distribution 
at angular distances larger than 0.15". The fitted 
gaussian is shown as the solid curve in Figure [3] 



A gaussian function was chosen since it provides a 
reasonable parametric representation of the dust 
emission. Interpolating this gaussian fit to the 
center of the disk suggests an expected central 
surface brightness of 31 mJy/beam, which is 4tj 
higher than the measured value. The significance 
level of the two intensity peaks, the fact that they 
appear in the map before cleaning, their orienta- 
tion along the disk major axis, and the symmetry 
with respect to the central star, suggest that they 
are real and, therefore, that the dust emission de- 
creases inside an orbital radius of about 14 AU. 
This is analogous to the situation in "transitional" 
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disks, where the inner gaps observed in the dust 
emission are att ributed to d usty dep leted inner re- 
ions (see, e.g.. lHughes et al. .2009i : iBrown et al. 



eions ( see, e.g. 
I2008il2009l) . 

At a first sight, this interpretation is incompat- 
ible with RY Tau's large near- and mid-infrared 
excesses, which suggest th e presence of warm dus t 
within 10 AU of the star (jRobitaille et al. 1120071 ). 
If, however, the inner disk is only partially de- 
pleted and dust emission remains optically thick 
in the infrared, the observed double intensity peak 
and the spectral energy distribution can be recon- 
ciled. A number of physical mechanisms could re- 
duce the dust density in the inner region of circum- 
stellar disks. For example, planets less massive 
than Jupiter may carve partially depleted gaps in 
the surface density distribution by tidal interac- 



tion w ith the surrounding material (jBrvden et al 
19991 ). This possibility is discussed in more de- 
tail in Section 15.31 In Paper I, we also proposed 



that a surface density profile that gradually de- 
creases towards the star may originate naturally 
from the viscous evolution of a disk if viscosity de- 
creases with radius. Finally, it is also possible that 
the decrease in dust emission may originate from 
a lowering of opacity du e to the growth of dust 



grain s to centimeter sizes (iDullemond &: Dominik 



20051 ). Unambiguously disentangling these models 



requires even higher angular resolution observa- 
tions than are yet available. 

We must note that radial velocities stud- 
ies ( Herbig fc Bell 19981) and Hipparcos obser- 
vations of the vari ability of the photocenter 
(|Bertout et al. I Il999l ) suggest that RY Tau is a 
binary. The Hipparcos data implies a minimum 
projected separation of 23.6 mas and a position 
angle of 304°±34°, almost perpendicular to the 
position angle of the disk inferred from our obser- 
vations. Assuming that the binary and the disk 
have the same inclination, the spatial separation 
between the binary components is 6-9 AU, and 
could explain the double peak in the dust contin- 
uum emission. Indeed, the presence of a stellar 
mass companion orbiting at a radius of 6-9 AU 
would push the inner radius of the circumstellar 
disk to a distance of 9-13 AU by tidal interactions 
( Wolf et al. II2007I) . However, the binary nature of 
RY Tau has been rendered questionable by near- 
infrared interferometric observations that suggest 
an inner disk radius at 0.1 AU from the central star 



and exclude the presence of a stellar mass compan- 
ion between 0.35 AU and 4 AU down to a stellar 
flux r atio of 0.05 (|Akeson et al . 2005; Pott et al. 
20091) . A stellar companion was also undetected 



in recent spectroscopic and aperture masking ob- 
servations (Duy Cuong Nguyen and Adam Kraus 
private communication). As discussed above, the 
spectral energy distribution is also inconsistent 
with the existence of a large inner gap completely 
depleted of gas and dust as could be expec ted for a 
stellar companion (i Robitaille et al. II2007I) . These 
results suggest that RY Tau is indeed a single star, 
and the variability observed by Hipparcos and the 
radial velocity variations may be attributed to 
brightness changes in the circumstellar environ- 
ment (see the discussion in Schegerer et al. 20081 
and references therein). 

A notable characteristic of our images of the 
dust emission is the high degree of central sym- 
metry and, with the exception of the innermost 
region, the almost complete absence of features in 
the surface brightness distribution. If the emission 
is optically thin (we will examine this assumption 
is Section [573]) . this translates to a smooth radial 
profile for the dust. The degree of symmetry of the 
emission can be quantified by analyzing the imag- 
inary part of the correlated fiux, plotted in Fig- 
ure [2] as a function of the angular frequency B^v ■ 
Point symmetric emission will have a zero imagi- 
nary part at all spatial frequencies. For RY Tau, 
the deviations from zero are comparable to the 
noise in the observations (see the left panels of 
Figure [2]) 

3.2. DG Tau disk morphology 

The surface brightness distribution for the 
DG Tau disk has a central peak coincident with 
the stellar position and smoothly decreases out- 
ward to reach the noise level at an angular dis- 
tance of about 0.5". At both 1.3 and 2.8 mm 
the emission appears fairly symmetric, and in- 
deed the imaginary part of the correlated flux is 
zero for Buv > 300 kA (see upper right panel of 
Figure [2]). The imaginary part increases up to 
50 mJy at shorter spatial frequencies, suggesting 
that the emission may be asymmetric on angular 
scales larger than l/B^v ~ 0.7". Comparing the 
real and imaginary parts of the correlated flux at 
the shortest spatial frequencies, we find that the 
asymmetric part of the emission contributes about 
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Fig. 3. — : Radial profile of the 1.3 mm surface 
brightness in RY Tau disk measured along the ma- 
jor axis of the disk at a position angle of 24° East 
from North. The intensity error bars (red points) 
correspond to the noise level of 0.9 mJy/beam. 
The solid curve shows a gaussian fit, while the 
dashed line indicates the zero intensity level. The 
bracket in the upper left indicates the angular res- 
olution. 



14% to the total flux. 

As already noted, the high angular resolution 
observations of RY Tau and DG Tau allow us to 
examine disk properties that bear on planet for- 
mation. In particular the radial density profile 
of the circumstellar material is critical to under- 
standing where planet formation may occur, or 
where it has occurred. With observations at more 
than one wavelength, we may also consider radial 
variations of the grain properties. As we describe 
below, measurements of the dust radial profile and 
the variation of the dust properties with radius 
are best undertaken by comparing the observa- 
tions with theoretical disk models in the Fourier 
domain. There, the effects of different angular res- 
olutions, sensitivity and atmospheric seeing can be 
more easily taken into account. Moreover, com- 
parison with theoretical models is the only way 
to quantify the contribution from optically thick 
emission, leading to an improved estimate of the 
dust mass. 

4. Disk and dust models 

To investigate the dust radial distribution 
around DG Tau and RY Tau, we consider two 



different models for the disk surface density. The 
first model consists of the classical power law pa- 
rameterization 



^40 



V40Auy 



for Rin < R < Ro 



(1) 

where S40 is the disk surface density at a ra- 
dius of 40 AU. Rin and Rout si's the inner and 
outer disk radii respectively. The second model is 
the similarity soluti on for the evolution of a vis- 
cous Keplerian disk (Lynden-Bell &: Pringle Ill974 : 



Hartmann et al. 11199 
this has the form 



As discussed in Paper I, 



X exp ■ 



1 



2(2-7) 



R 



(2-7) 



(2) 



where is the surface density at radius Rt , some- 
times called the transition radius. For R < Rt, the 
surface density has a power law profile character- 
ized by the slope 7, while at larger radii the surface 
density falls exponentially. 

These two different parameterizations are used 
to calculate the dust emission by solving the struc- 
ture of an hydrostatic disk heate d by the stellar 
radiation ( Dullemond et al71l200l[ ). The tempera- 
ture on the disk mid plane is self-consistently cal- 
cu lated by adopting the two la yer approximation 



of IChiang fc Goldreich I (Il997h . The disk tem- 



perature, which depends mainly on the disk op- 
tical depth at optical and infrared wavelengths, 
is roughly proportional to for both surface 

density profiles (see Paper I and reference therein 
for a detailed discussion on the disk temperature 
radial profile). 

Fundamental to any disk model is the prescrip- 
tion adopted for the dust opacity. Although the 
dust contributes only about 1% to the total disk 
mass, it dominates the disk opacity in the wave- 
length range where most of the stellar and disk ra- 
diation is emitted. We assume that the dust size 
distribution follows a power law n{a) cx a~' where 
a is the radius of a dust grain. The assumptions on 
the slope q, on the minimum and maximum grain 
sizes, on the dust chemical composition and on the 
grain structure define the frequency dependence 
of the dust opacity /c^, and, ultimately, the disk 
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Table 2 

Properties of the adopted dust models 



Object O^ynax Q ^1.3 mm 

(cm) (cm^/g) 


mm 

(cmVg) 




High dust opacity model, H 
DG Tau 0.075 3 0.082 
RYTau 0.035 3 0.131 


0.056 
0.075 


Low dust opacity model, L 
DG Tau 5.0 3.5 0.012 
RYTau 5.0 3.7 0.026 


0.0078 
0.015 



emission. The dust opacity is calculated for com- 
pact spheres composed o f astronomical silicates 



and organic carbonates (jWeingartner fc Draine 



200 It IZubko et al. 1 119961 ) . We assume a mass ra- 



tio of 1 between silicates and organics, which leads 
to grain density of 2.5 g/cm^. The dust opacity 
averaged over the grain size distribution is calcu- 
lated by fixing the minimum grain size to 0.005 
lim. The maximum grain size amax and the slope 
q are set to reproduce the observed slope of the 
spectral energy distribution as discussed below. 

At millimeter wavelengths the dust opacity can 
be approximated by a powe r law kn = A:o(i^/t'o)^ 
( Beckwith fc Sargent 199l[ ). If the dust emission 
is optically thin and the Rayleigh- Jeans approxi- 
mation is satisfied, the slope /3 of the dust opac- 
ity is related to the spectral index a of the ob- 
served disk emission F^, {F^, oc p") by the relation 
a — 2 + (3. This relation is only approximate if the 
dust emission is optically thick at some radii. In 
Paper I, we derived values for /3 of 0.5 and 0.7 for 
DG Tau and RY Tau respectively from an analysis 
of the SED, taking into account the optically thick 
contribution to the total dust emission. For the as- 
sumed dust composition and structure, these val- 
ues of /3 can be reproduced with different choices 
of the maximum grain size a^ax smd the grain 
size slope q (see Appendix E| . To investigate how 
the assumptions on the grain size distribution af- 
fect the model fitting, we adopt two different dust 
models that correspond to the extreme cases of 
low (L) and high (H) opacity. The correspond- 
ing dust opacities at both 1.3 mm and 2.8 mm are 
given in Tabled 

Finally, we assume that the dust opacity is con- 



stant throughout the disk. This is indeed one of 
the main assumption we want to test by modeling 
the observed dust emission at 1.3 mm and 2.8 mm 
and will be discussed in detail in Section [5^ 



5. Results and Discussion 

Models and observations are compared in 
Fourier space to avoid the non linear effects in- 
troduced by the cleaning process. The best fit 
models are found by minimization with five 
free parameters: the disk inclination i, the disk 
position angle PA, Rout, S40, and p for the power 
law surface density (Equation [1]), and i, PA, Rt, 
St and 7 for the similarity solution (Equation [2|) • 
The disk inner radius Rin is fixed at 0.1 AU. For 
both surface density models we find best fit so- 
lutions for both the high (H) and low (L) dust 
opacity models. The 1.3 mm and 2.8 mm data are 
fitted independently. 

To minimize ^ind evaluate the constraints 
on the model parameters, we use a Bayesian ap- 
proach that adopts uniform prior probability dis- 
tributions. In practice we sample the proba- 
bility distribution by varying the free parameters 
using the Markov Chain Monte Carlo method de- 
scribed in Paper I. 

Once a best fit solution is found, we confirm 
that this indeed corresponds to an absolute mini- 
mum of x^, as opposed to a local minimum, by 
running multiple Monte Carlo simulations with 
random initializations and verifying that they all 
converge to the same solution. Each parame- 
ter is allowed to vary in a large range: 0-80° 



8 



Table 3 

Best fit parameters assuming H dust model. 



Similarity solution 



Object 


A (mm) 


^n 


PA (°) 


Rt (AU) 


7 


St (g/cm^) 




DG Tau 


1.3 


24±9 


119±23 


23.4±1.8 


0.33±0.15 


10.9±1.5 


1.0608 




2.8 


31±12 


144±19 


27.7±3.0 


0.10±0.24 


7.5±1.3 


1.0629 


RYTau 


1.3 


66±2 


24±3 


26.7±1.2 


-0.54±0.18 


2.6±0.2 


1.0896 




2.8 


71±6 


20±4 


26.5±2.7 


-0.08±0.54 


2.6±0.5 


1.1894 


Power law 


Object 


A (mm) 




PA (°) 


Rout (AU) 


P 


i;4o (g/cm2) 




DG Tau 


1.3 


27±8 


120±26 


72.6±6.3 


1.00±0.15 


5.6±1.5 


1.0611 




2.8 


32±11 


144±18 


82.2±10.5 


0.74±0.24 


4.5±1.6 


1.0629 


RYTau 


1.3 


66±2 


24±3 


70.6±3.9 


0.12±0.15 


1.9±0.6 


1.0897 




2.8 


71±6 


20±4 


76.9±12.0 


0.64±0.45 


1.6±1.0 


1.1894 



Note. — The uncertainties correspond to a likelihood of 99.7% (i.e. 3a) for the normal distributions shown 
in Figure [SHll 



for the inclination, ±90° for the position angle, 
10 - 1000 AU for Rt and Rout, ±4 for p and 7, 
and 0.1 - 1000 g/cm^ for S40 and E*. 

The best fit disk models found for high and low 
dust opacities are listed in Table [3] and Table H] re- 
spectively. Each table lists the parameters for the 
similarity solution disk model in the upper part, 
and for the power law disk model in the lower part. 
The probability distributions for each free param- 
eter are shown in Figure|4]and Figure[5]for RY Tau 
in the case of the similarity solution and power law 
respectively. The same quantities for DG Tau are 
shown in Figure |6] and Figure [T) In these figures, 
the black and red histograms indicate the proba- 
bility distributions derived by fitting the 1.3 mm 
and 2.8 mm observations, respectively; solid and 
dashed curves represent the H and L dust opacity 
models. For each parameter we derive the uncer- 
tainty range that corresponds to a likelihood of 
99.7% (3(t) by fitting a normal distribution to the 
probabilities. 

Finally, Figure [5] shows comparisons between 
the observed real part of the correlated flux (filled 
squares with error bars), the best fit models for 



the similarity solution (solid curve), and a power 
law surface density (dashed curve). 

5.1. Dependence on the dust opacity and 
implications on the disk masses 

The best fit solutions for the H and L dust 
opacity models are shown in Figure (SUB] with solid 
and dashed curves respectively. In all cases, H 
and L models lead to very similar values for the 
disk position angle, the disk inclination and the 
radial profiles of the surface density defined by 
p and Rout in the case of the power law mod- 
els, and 7 and Rt for the similarity solution mod- 
els. As discussed in Paper I, these parameters 
are essentially independent of the dust opacity. 
This is mainly because the disk mid-plane tem- 
perature Ti(R) varies by only a few percent be- 
tween the different dust models, as long as the disk 
is optically thick to the stellar radiation. Since 
I](i?) cx Ti{R)^^, the radial profile of the surface 
density varies by only small fraction when different 
dust models are assumed. 

By contrast, the surface density normalization 
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Table 4 

Best fit parameters assuming L dust model. 



Similarity solution 



Object 


A (mm) 




PA (°) 


Rt (AU) 


7 


St (g/cm^) 




DG Tau 


1.3 


24±11 


119±24 


22.5±1.8 


0.28±0.15 


74.4±9.9 


1.0608 




2.8 


31±12 


144±20 


26.4±2.7 


0.07±0.27 


55.4±8.9 


1.0629 


RYTau 


1.3 


66±2 


24±3 


25.6±1.2 


-0.58±0.18 


13.6±1.2 


1.0896 




2.8 


71±6 


20±4 


25.1±2.4 


-0.10±0.57 


14.3±2.3 


1.1893 


Power law 


Object 


A (mm) 


^n 


PA (°) 


Rout (AU) 


P 


i;4o (g/cm2) 




DG Tau 


1.3 


27±9 


120±24 


72.3±4.0 


1.06±0.18 


35.7±3.6 


1.0611 




2.8 


32±11 


144±19 


81.8±9.3 


0.74±0.24 


32.1±4.5 


1.0629 


RYTau 


1.3 


66±2 


24±3 


70.5±3.9 


0.11±0.18 


9.7±1.2 


1.0897 




2.8 


71±5 


20±4 


76.7±12.6 


0.68±0.51 


8.3±2.3 


1.1893 



Note. — The uncertainties correspond to a likelihood of 99.7% (i.e. Za) for the normal distributions shown 
in Figure [SHll 



(E( and S40) varies with the dust opacity so that 
the product S x /c^ remains almost constant if the 
emission is optically thin. Consequently, a lower 
dust opacity requires a higher dust mass in order 
to emit the same amount of radiation at millime- 
ter wavelengths. The ratio E^/Sj^ is then ap- 
proximately equal to the ratio between the dust 
opacities listed in Table [H 

From the analysis of the surface density of the 
best fit model we find that the RY Tau emis- 
sion is always optically thin at both 1.3 and 
2.8 mm. However, DG Tau emission is optically 
thick within 20 AU at 1.3 mm for both the sim- 
ilarity solution and the power law models. The 
1.3 mm flux emitted within this region is about 
25% of the total flux. At 2.8 mm the emission is 
always optically thin in the case of the similarity 
solution while it is optically thick within 6 AU in 
the power law case. In this case the optically thick 
contribution is 5% of the total flux. 

Different dust opacities lead to different values 
for the total mass of dust in the disks. For DG Tau 
we obtain total dust masses of about 33 and 233 
Earth masses (A/^) in the case of the high and 



low opacity dust models respectively. Disk masses 
of ~10 and 50 are found for RY Tau. Mas- 
sive disks can also be obtained by extending the 
grain size distribution larger than 5 cm. For ex- 
ample, in Paper I we derived total dust masses of 
about 1331 and 216 for DG Tau and RY Tau 
respectively by assuming a maximum grain size 
of 10 cm and a slightly different grain composi- 
tion. Additional uncertainties in the disk mass 
come from the dust chemical composition. As dis- 
cussed in Appendix El the presence of ice or vac- 
uum in the grains leads to smaller dust opacities 
at millimeter wavelengths and consequently pro- 
duces higher disk masses. We therefore estimate 
that the circumstellar disks around DG Tau and 
RY Tau contain a minimum mass of dust of 30 and 
10 M0 respectively, while the upper limit is not 
constrained due to the uncertainties on the grain 
size distribution. Assuming the standard dust/gas 
ratio of 0.01, these values correspond to minimum 
disk masses of 0.009 and 0.003 Mq for the DG Tau 
and RY Tau respectively. 
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Fig. 4. — : Probability distribution of the disk pararaeters obtained by fitting the RY Tau observations at 
1.3 mm and 2.8 mm with the similarity solution for the surface density distribution. The results obtained 
by fitting the 1.3 mm data only are shown by the black curves, while the red curves indicate the results 
obtained by fitting the 2.8 mm only. Solid curves show the probability distribution obtained assuming the 
dust opacity model H, and dashed curves correspond to the dust opacity model L (see Table [5]). 



5.2. Constraints on the surface density: 
similarity solution vs power law 

Figure |8] shows the comparison between models 
and observations in terms of the real part of the 
correlated flux as a function of the baseline length. 
To correct for the disk inclination we deprojected 
the baseline assuming the inclinations and position 
angles listed in Table [31 In this figure, the results 



for H and L dust models lead to indistinguishable 
curves. Similarity solution and power law mod- 
els are represented with solid and dashed curves 
respectively, and the observations are shown by 
black dots with error bars. It is clear that both 
the similarity solution and power law disk models 
provide satisfactory fits to the observations. The 
similarity solution model provides smaller values 
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Fig. 5. — : Probability distribution of the disk parameters obtained by fitting the RY Tau observations at 
1.3 mm and 2.8 mm with a power law surface density distribution. Colors, solid curves, and dashed curves 
are the same as in Figure HI 



of (see Table [3] and S]) and, in the case of DG 
Tau, a better fit to the observations between 400 
and 800 kA. In this range of spatial frequencies, 
the power law solution is characterized by a wiggle 
due to the sharp truncation of the dust emission 
at 72 AU. On the other hand, the exponential ta- 
pering of the similarity solution leads to a smooth 
visibility profile that matches extremely well the 
observations. The same behavior is present in the 
lower panel which compares the model and the 



observations at 2.8 mm. However, in this case the 
observations at Buv > 400kA are too to distinguish 
between the two models. Although not conclusive, 
this result make the similarity solution model a 
more appealing explanation for the dust emission 
in circumstellar disks, confirming the conclusions 
of iHughes etai] (|2008l ). 

Figure ini shows the surface density derived from 
the 1.3 mm observations for both the power law 
and the similarity solution model in the case of 
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Fig. 6. — : Probability distribution of the disk parameters obtained by fitting the DG Tan observations at 
1.3 mm and 2.8 mm with the similarity solution for the surface density distribution. Colors, solid curves, 
and dashed curves are the same as in Figure 



high dust opacity. The two models lead to sim- 
ilar values of S(-R) in the region where most of 
the 1.3 mm flux is emitted, namely between ~15 
and 50 AU. In this region, the surface density in 
RY Tau disk is almost constant with the radius, 
while it decreases roughly as 1/i? in the case of 
DG Tau. Inside 15 AU and outside 50 AU, the 
observations lack both the angular resolution and 
the sensitivity required to directly constrain the 
surface density. As a consequence, the values of 



strongly depend on the assumed model and 
can differ by one order of magnitude at the disk 
inner radius. 

5.3. Surface density and implication on 
the existence of planets 

In this section we discuss the implications of the 
inferred surface density on the presence of planets. 
The analysis is limited to surface density profiles 
obtained by fitting the observations at 1.3 mm. 
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Fig. 7. — : Probability distribution ol the disk parameters obtained by fitting the DG Tau observations at 
1.3 mm and 2.8 mm with a power law surface density distribution. Colors, solid curves, and dashed curves 
are the same as in Figure HI 



which have the highest angular resolution. 

5.3.1. DG Tau 

For the similarity solution model, the surface 
density has a radial profile characterized by 7 ~ 
0.31 ± 0.18 and i?t - 23 ± 2 AU. The transition 
radius Rt agrees well with our earlier observations 
(21±3) but 7 is significantly larger than the value 
of -0.5±0.6 from Paper I. The discrepancy is prob- 
ably due to the fact that the earlier observations 



were taken in poorer weather conditions and the 
model fitting did not account for the atmospheric 
seeing. Figure 10 shows the residuals after sub- 
tracting the best fit model to the new observations. 
Note that the power law model gives very similar 
residuals. The residuals are as high as 3-6cr and 
are found at angular scales larger than 0.7" where 
the emission is slightly asymmetric (see Section[3|). 
In this outermost disk region, the surface density 
may deviate significantly from the symmetric ra- 
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Fig. 8. — : The real part of the correlated flux as function of the deprojected baseline. The solid curves 
show the best-fit models where the surface density is described by the similarity solution for viscous disk 
evolution. The dashed curves show the results for a model with a power law surface density. The models 
were computed using the high dust opacity, but identical results were obtained for the low dust opacity case. 



dial profile assumed in the model. We calculate 
that variations of ±10 — 30 g/cm^ with respect to 
the best fit surface density profile over a spatial 
region comparable with the beam size may pro- 
duce the observed residuals. Larger variations of 
the surface density on smaller angular scales are 
also possible. 

The residuals do not show global deviations 
from the smooth surface density profile, appar- 
ently excluding the possibility of gaps in the dust 
distribution that might be produced by a planet. 
Of course, low mass planets may not produce any 
discernible gap and may still exist in the DG Tau 
disk. The formation of a gap is possible only if 
the efficiency in removing the material close to 
the planet orbital radius via tidal torques is larger 



than the mass a ccretion rate due to the disk vis- 
cosity (see, e.g, iLin fc Papaloizou 1 119931 ) . If we 



assume the a prescription for the disk viscosity 
( Shakura &: Sunvaev 1 19731) . and call h the pres- 
sure scale height of the disk, a planet orbiting 
at radius Rp can open a radial gap in the disk 
surface density on l y if Mp/M^ > 32a; ( h/Rr,) ^ 
(|Lin fc Papaloizou I Il993t fervden et all Il999l ). 
Moreover, for a disk in hydrostatic equilib- 
rium with the gravitational field of the cen- 
tral star, the pressure scale h is proportional 



to 



T,iRp)^/^R;'^M: 



■3/2s,f-l/2 



(|Chiang fc Goldreich 



1997h . Since the temperature is Ti{Rp) oc Rp 
(see Paper I for more details), the formation of a 
gap requires 



Mp > 26.3i?y2a, 



(3) 
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Fig. 9. — : Surface density for DG Tau (top panel) 
and RY Tau (lower) disk. The red and blue curves 
show the best fit solutions for the power law and 
the similarity solution models respectively in the 
low opacity case L. The dashed curves show the 
3cr uncertainty range for the surface density. The 
surface densities for the high opacity case will have 
similar shapes, but will be about a factor of 6 
lower. 



where Mp is expressed in Jupiter masses, Rp is 
AU, and the numerical constant is calculated for 
a disk temperature of 194 K at 1 AU as determined 
from our disk model. Typical values of a are in 
the range 10^^ — 10""^, and imply that a planet 
can open a gap at 1 AU only if its mass is larger 
than about 0.1 Mj. To open a gap at 30 AU, the 
mass must be larger than about 0.5 AIj. 

To investigate the effects that a planet more 
massive than 0.1 Mj might have on the obser- 
vations of the dust continuum emission, we sim- 
ulated the presence of a planet in the DG Tau 
disk by opening a gap in the surface density dis- 
tribution corresponding to the best fit models dis- 
cussed above. For simplicity, we assumed that 
the planet describes a circular orbit and that the 
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Fig. 10. — : The black curves show the residuals 
for the 1.3 mm DG Tau observations after sub- 
tracting the best-fit model for the similarity solu- 
tion. Contours start at 3cr and are spaced by la. 
The thin red curves and the color scale show the 
observed dust emission, with contours spaced by 
3cr. 



gap can be represented by a circular ring. To be 
compatible with numerical simulations of planet- 
disk interaction, the half-width of the ring A 
is assumed to be equal to t wice the Hill radius 

119991: 



Rh = Rr,VMp/ {2,M^,) (e.g. iBrvden et al 



Wolf et al.l[2007h . In the region between i?p ± A 



the surface density is depleted by a fraction / 
that depends on the mass of the planet and on 
the disk viscosity. For a = 10"'^, we can assume 
/ = for planet masses Mp > \Mj, f = 0.1 



/ = 0.17 for M, 



0.3 M 



0.1 Mj (IWolf et al. 1120071) 



J 



for Mp = 0.5 Mj, 
and / = 0.6 for Mp 
Therefore, only planets more massive than 1 Mj 
will produce completely cleaned gaps. 

We simulated gaps corresponding to planets in 
the mass range 0.3-5 Mj and with orbital radii 
between 1 and 90 AU. For each model we calcu- 
lated the residuals as the difference between the 
observations of DG Tau at 1.3 mm and the model 
image. If the gap is too small compared to our 
angular resolution, or too faint compared with our 
sensitivity, the residuals will be similar to the case 
without gaps shown in Figure 1101 In this case we 
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say that the gap is not detected. On the other 
hand, a large and deep gap will produce a bright 
ring in the residual. To quantify how reliable the 
detection of a gap is, we define a signal-to- noise 
ratio of the gap, gap SNR, in the following way. 
First, we deproject the residual for the inclination 
and position angle of DG Tau disk. Then we take 
the radial average of the residuals at the distance 
corresponding to the orbital radius of the planet 
adopting a radial bin width equal to the FWHM of 
the synthesized beam (i.e., 0.17"). We define the 
gap SNR as the mean residual in the radial bin 
divided by the uncertainty in the mean. In this 
way, detected gaps correspond to gap SNR> 3. 

The results are summarized in the upper panel 
of Figure [TTJ Planets with masses and radii that 
lead to gap SNR> 3 produce detectable gaps. 
No gaps are detected in our observations of the 
DG Tau disk at more than 3cr (see the red curve). 
This enable us to constrain the masses and orbital 
radii of any planets that may be present. In par- 
ticular, we can exclude that planets more massive 
than Jupiter exist between 5 and 40 AU, or that 
planets with masses slightly smaller than Jupiter 
exists between 10 and 25 AU. The observations 
lack both the angular resolution and the sensitiv- 
ity required to detect gaps produced by planets 
with a mass smaller than about 0.5 Mj. 

An important caveat is that a planet may exist 
but may not have had enough time to completely 
open a gap in the disk. The gap formation time 
scale TA results from the trade-off between the ef- 
ficiency of the tidal torque exercised by the planet 
in removing angular momentum, and the accre- 
tion of new material coming from larger radii in 
the gap due to the disk viscosity. A lower limit 
of the gap formation time scale is obtained in the 
zero viscosity limit. In this ca se, an an a lytic for- 



mulation is provided by Brvden et al. 
the form 

P / A 



'A 



g2 Vi?; 



(4) 



where P is the orbital period, q — Mp/Mi, and 
A — 2Rh as defined above. Assuming Keplerian 
rotation, we can rewrite the time scale for the gap 
formation as 



l.lMyrx 



3/2 



Rp 
AU 



3/2 



Mj 



A 

R„ 



(5) 




30 40 50 60 
Orbital radius (AU) 



90 



Z 

o 




10 20 30 40 50 60 
Orbital radius (AU) 



90 



Fig. 11. — : The signal-to-noise ratio of the detec- 
tion of a gap generated by a planet as a function 
of the orbital radius Rp and the planet mass. The 
different curves correspond to masses between 0.5 
and 10 Jupiter masses as labeled. The thick red 
curves indicate the signal-to-noise ratio measured 
from the 1.3 mm images after subtracting the best 
fit model for the similarity solution. 

For Gap SNR> 3, planets should produce a 
detectable gap. The upper and lower panel refers 
to the case of DG Tau and RY Tau respectively. 



The upper panel of Figure [12] shows the calcu- 
lated values of t^*" for the stellar mass of DG Tau 
(0.3 -M0). In the case of a planet with a mass be- 
tween 0.3 and 0.5 Mj orbiting at a radius larger 
than 40 AU, the minimum time scale for the gap 
formation is comparable with the age of the system 
(0.1 Myr). For more massive planets, or for closer 
radii, the minimum gaps time scale is a small frac- 
tion of the age of the system. 

We conclude that, for DG Tau, the observations 
lack the sensitivity and angular resolution required 
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Fig. 12. — : Minimum time scale (rA)for the for- 
mation of a gap as a function of the orbital radius 
and mass of the planet. The different curves cor- 
respond to planet masses between 0.3 and 5 Mj as 
labeled in the figure. The upper and lower panels 
show the time scale for a central star of 0.3 Mq 
and 2 Mq respectively. 

to investigate the presence of planets less massive 
than about 0.5 Mj at any orbital radius. Our 
analysis indicates that no planets more massive 
than Jupiter are present between 5 and 50 AU, 
unless they are younger than 10^ years. 

5.3.2. RY Tau 

The similarity solution for the disk surface den- 
sity is characterized by 7 = —0.56 ± 0.18 and 
i?( ~ 26±3 AU. As shown in FigurelHl the surface 
density increases roughly as ^/R from the inner 
radius at 0.1 AU up to about 26 AU and then 
decreases exponentially outward. This supports 
the suggestion in Section [3] that the RY Tau in- 
ner disk might be partially dust depleted with re- 



spect to power law disk models. We note that 
this surface density profile may provide an expla- 
nation for both the double peak intensity at 1.3 
mm and the disk excess at infrared wavelengths. 
Indeed, within 10 AU the model disk remains op- 
tically thick at optical and infrared wavelengths, 
exhibiting the infrared excess typical of classical 
disks. 

At larger radii, the surface density in RY Tau 
disk decreases smoothly and the residuals calcu- 
lated by subtracting the best fit models to the 
1.3 mm dust emission map do not show any struc- 
ture at more than 3(t. This excludes strong devi- 
ations from an unperturbed viscous disk profile. 

The lower panel of Figure [TT] shows the signal- 
to-noise ratio of the detection of a gap generated 
by planets of 1, 5 and 10 Jupiter masses as a 
function of the orbital radius. Due to the higher 
disk inclination and stellar mass, a planet orbit- 
ing around RY Tau would produce a less visible 
gap. In particular, our observations seems to ex- 
clude the presence of planets more massive than 5 
Jupiter masses between 10 and 60 AU. Given the 
higher stellar mass of RY Tau, the minimum time 
scale for the formation of gaps is one order of mag- 
nitude larger than the case of DG Tau (see lower 
panel of Figure [T2)) . This implies that planets less 
massive than Jupiter orbiting at more than about 
30 AU may not have had enough time to form a 
gap in the disk. 

5.4. Radial dependence of the dust prop- 
erties 

A comparison of the best fit solutions obtained 
for the wavelengths of 1.3 mm and 2.8 mm enable 
us to investigate the dependence of dust opacity 
on the orbital radius. If the dust opacity is con- 
stant throughout the disk as assumed in SectionHJ 
the model fitting necessarily leads to the same sur- 
face density profile for observations at two differ- 
ent wavelengths. Otherwise, different would 
suggest a radial variation in the relative dust opac- 
ities at the observed wavelengths. To understand 
this point, we assume that the dust emission is op- 
tically thin. In this case the observations constrain 
the product Y,x{R) x fc^, where Yix{R) is the sur- 
face density obtained by fitting the observations 
at the wavelength A. In the more general case in 
which the dust opacity depends on the orbital ra- 
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dius we can write 



(6) 



The left side of this equation contains the opac- 
ity discussed in Section |4] and the surface density 
derived from the model fitting. The right side con- 
tains the unknown "true" surface density E(i?) in 
the case in which the "true" dust opacity k\{R) 
varies with the radius. The ratio of Equation |6] for 
two different wavelengths Aq and Ai leads to 



/3(fl) 



(7) 



Here we assumed that at each radius the dust 
opacity can be expressed by a power law k\ oc . 
Finally, taking the logarithm of this latter equa- 
tion we can write 



/3{R) ^(3, + Af3{R) 



(8) 



where /3c = log(fcAi //cao)/ log(Ao/Ai) and A/3(i?) 
has the form 



A/3(i?) = log- 



X log 



^xAR) 



(9) 



If Eai = 5^Ao I the dust opacity slope is constant 
throughout the disk and assumes the value dis- 
cussed in Section ID Otherwise, we can use the 
latter equation to investigate the radial variation 
of f3. 

The best fit solutions for obtained at 

1.3 mm and 2.8 mm are shown in Figure [5][6] with 
black and red curves respectively. The best fit pa- 
rameters are summarized in Table [3] and IH The 
quoted uncertainties correspond to a likelihood of 
99.7% (i.e. 3a) and are calculated by fitting a nor- 
mal distribution to the measured probability dis- 
tributions. For RY Tau, the disk model obtained 
by fitting the two wavelengths separately are in 
agreement within Scr. For DG Tau the solutions 
disagree by more than 3ct only in the case of the 
similarity solution and high dust opacity. 

Figure [T3l shows the radial variation of /3 as de- 
fined in Equation [8] for both DG Tau and RY Tau. 
The region marked with color indicates values of 
j3 within 3ct from the radial profile corresponding 
to the best fit solution for the surface density in 
the case of the similarity solution model. Values 
of /3 outside this region are rejected by our obser- 
vations. 



The results for RY Tau and DG Tau are clearly 
consistent with a large variety of radial profiles of 
/?. For both sources, /? is better constrained be- 
tween radii of 20 and 70 AU. However, even in this 
interval, the observations constrain possible vari- 
ation of /3 to within only A/3 < 0.7. Nevertheless, 
across most of the disk the circumstellar dust dif- 
fers from that observed in the interstellar medium 
(ISM). Dust in the ISM is characterized by sub- 
micron dust grains and by a millimeter opacity 
slope of /3 ~ 1.7. In contrast, in both sources, /3 
is smaller than 1.7 up to at least 80 AU, suggest- 
ing that the circumstellar dust has been processed 
and, in particular, has increased its size up to a 
maximum value that varies between 20 ^m and a 
few centimeters. 

Although in both sources /3 may be constant 
throughout the disk (see dashed lines), our results 
suggest that /3 decreases with the radius in DG 
Tau disk. As discussed in Section |31 /3 depends on 
a number of poorly constrained quantities, such as 
composition, structure, and size of the dust grains. 
For example, varying /3 from ~1 to ^0.2, similar 
to what is suggested for DG Tau between 10 and 
60 AU, may be due to the maximum grain size 
increasing from 20 fim to 1 cm for (7 = 3, or, alter- 
natively, to a decrease of q from 4 to 3 if the max- 
imum grain size is between 1 and 10 cm. In short, 
the interpretation of /3 only in terms of the grain 
size distribution can be very misleading. It seems 
most plausible that both dust composition and the 
relative contributions of smaller and larger grains 
change through the disk, contributing to the vari- 
ation of the dust opacity. It is clear that better 
constraints on the radial profile of /3 are required 
before pushing the investigation further. 

In this regard, we note that the current con- 
straints on the radial variation of /3 are limited by 
two factors. First, although the angular resolution 
of the observations described here is significantly 
better than hitherto possible, the dust surface den- 
sity is well constrained only between 15 and 50 
AU, where most of the observed flux is emitted 
(see the discussion in Section At smaller and 
larger radii the surface density is uncertain by al- 
most one order of magnitude. Second, our analysis 
is hampered by the small separation in wavelength 
between the observations since A(3{R) is propor- 
tional to log~^{Xi/ Xf)). The uncertainties shown 
in Figure [T^] can be reduced by a factor of 2 by 
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Fig. 13. — : Slope of the millimeter dust opacity /3 as a function of the radius for RY Tau (left panel) and 
DG Tau (right panel). Values of f3 outside the colored region are excluded by our observations at more 
than the 3cr level. The dashed lines show the mean value of /3 derived in Paper I as discussed in Section |4l 
The thick solid curve shows the radial profile of /3 corresponding to the best fit solution for the disk surface 
density at both 1.3 and 2.8 mm. The solid horizontal line shows the value of /3 characteristic of the dust in 
the interstellar medium. 



extending the observations at 7 mm. These obser- 
vations will become possible with the expanded 
correlator on the EVLA. 

6. Conclusions 

We have presented CARMA observations of 
the dust thermal emission at the wavelengths of 
1.3 mm and 2.8 mm from the circumstellar disks 
around the pre-main sequence stars RY Tau and 
DG Tau. The observations are characterized by 
unprecedented angular resolution of ~0.15" and 
0.30" at 1.3 mm and 2.8 mm respectively, cor- 
responding to spatial scales of 20 and 40 AU at 
the distance of Taurus. Based on these images, 
we have addressed three fundamental questions re- 
lated to the formation of planets in the disk around 
pre-main sequence stars. What is the radial den- 
sity distribution of circumstellar dust? Does the 
dust emission show any indication of the presence 
of (proto)-planets? Do the dust properties vary 
with orbital radius? 

By analyzing the morphology of the surface 
brightness of the dust emission and comparing 
the observations with theoretical disk models, we 
make the following conclusions: 

(1) Both the classical power law disk surface 
density (Hayashi 1981) and the similarity so- 



lution for the viscous evolution of a Keple- 
rian disk (Hartmann 1998) fit the observa- 
tions well. The surface density is well con- 
strained between 15 and 50 AU. In this re- 
gion, the two models lead to values of S that 
agree within 30% for a fixed dust opacity. At 
smaller and larger radii, the surface density 
depends on the assumed model and varies by 
almost an order of magnitude. We have ver- 
ified that the assumptions on the dust opac- 
ity have a small effect on the model fitting 
and, therefore, on the radial profile of the 
dust density. However, the total disk mass 
may vary by almost two order of magnitude 
for different dust compositions and grain size 
distributions. 

(2) The dust emission in DG Tau is mostly radi- 
ally symmetric. It is characterized by a sin- 
gle, central peak and smoothly decreases up 
to an angular distance of about 0.5". Theo- 
retical disk models reproduce the observa- 
tion very well, with randomly distributed 
residuals between 3 and 6a. No system- 
atic deviation from the similarity solution for 
the surface density of a viscous disk are ob- 
served. By simulating the presence of plan- 
ets in the disk via the gap in the surface 
density produced by tidal torques, we find 
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that the observations exclude the presence of 
planets more massive than Jupiter orbiting 
between 5 and 40 AU from the central star, 
unless the planets are very young (< 10'* yr) 
and have not had the time to open a gap 
in the disk. The observations lack both the 
angular resolution and sensitivity to inves- 
tigate the presence of planets less massive 
than about 0.5 Jupiter masses. 
For RY Tau, the dust emission is charac- 
terized by two peaks separated by about 
28 AU that suggest a decrease in the sur- 
face density, or dust opacity, within 14 AU 
of the central star. We found that the sim- 
ilarity solution for the disk surface density 
is characterized by a negative value of 7, 
and provides a reasonable explanation for 
the double peak intensity observed at 1.3 
m m. Depletion of milhmeter dust grains 
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values of the disk viscosity, or the presence 
of planetesimals, may produce the observed 
dust morphology. At larger radii, the dust 
emission shows a very smooth profile with 
no asymmetries or gaps. The lack of gaps in 
the disk suggests that any planets between 
10 and 50 AU are less massive than about 5 
Jupiter masses, or, as for DG Tau, are very 
young. 

(3) The best-fit models to the 1.3 mm and 
2.8 mm data were compared to investigate 
the radial dependence of the slope opacity /3, 
assuming that the dust opacity at millime- 
ter wavelengths is expressed by a power law 
A-*^. We can exclude cases in which /3 
varies by more than 0.7 within 70 AU. Nev- 
ertheless, between 20 and 70 AU, the disks 
around DG Tau and RY Tau are character- 
ized by values of j3 that are smaller than 
that found in the ISM. This implies that the 
dust has been reprocessed and has grown in 
size up to a radius of at least 20 /.tm. The 
investigation of the radial variation of j3 is 
still limited by the angular resolution and by 
the small separation in wavelength between 
the observations. In the future, ALMA and 
the EVLA will play crucial roles in the in- 
vestigation of the radial dependence of the 
dust properties by increasing the angular 
resolution and the interval in wavelength. 
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A. EfTects of the assumptions on the dust composition on the dust opacity and size 



The slope of the dust opacity /3 measured at millimet er wavelengths has been widely adopted to constrain 



the size and opacity of the circumstellar dust (see, e.g., Ricci et al. 2009f) . In this appendix we investigate 



how these latter two quantities depend on the assumption on the grain size distribution and composition. 

The color scale in Figure [Ml shows the dust opacity at 1.3 mm calculated for three different dust com- 
positions. In model A we assume that the grains are compact spheres c omposed of astronomical silicates 



and organic carbonates (jWeingartner fc Draine II2OOII : IZubko et al. 1119961 ). This is the dust model adopted 



in the paper. In model B the grain is composed of silicates, carbonates and water ice. Finally, in model C 
the ice is replaced by a vacuum, resulting in a porous grain made of silicates and carbonates. We assume a 
mass ratio of 1 betw een silicates and organics, and of 0.7 between silicates and the water ice or the vacuum 
( Pollack et al. The resulting bulk densities are 2.5 g/cm^ for model A, 1.9 g/cm'^ for model B and 

1.5 g/cm'^ for model C. We assume a grain size distribution n{a) oc a^' and fix the minimum grain size to 
0.005 yum. We then calculate the dust opacity by varying the maximum grain size Umax between 6 x 10~^ 
and 10 cm and the slope q between 2 and 5. 

The solid thick curves show the values of {amax,Q) required to obtain the values of /? measured for DG 
Tau and RY Tau, 0.5 and 0.7 respectively. From the figure, it is clear that the value of sets a lower bound 



Natta et al~ll2004[ ). This lower bound 



to the maximum grain radius in the grain size distribution (see, e.g. 
is a function of the grain composition and increases by almost an order of magnitude between model A and 
C. For example, in order to have /3 = 0.7, the maximum grain size must be at least 0.03 cm for compact 
grains in model A, or about 0.2 cm for the porous grains in model C. Note, however, that we can obtain the 
same value of /3 with amax — 10 cm and q between 3.2 and 3.7. Deriving the maximum grain size from the 
measure of (3 is evidently strongly degenerate. 

This introduces a large uncertainty on the dust opacity and ultimately on the total mass of circumstellar 
dust derived fro m millimeter observa tions. Even if we limit the analysis to the generally adopted value 
3.5 (see, e.g.. iBrauer et al. II2008I ). the dust opacity for (3 = 0.7 varies from about 8.4 cm^/g of dust in 
the case of the model A and Umax = 0.1 cm, to 1.2 cm^/g of dust for model C and Umax = 10 cm, leading 
to disk masses that differ by almost one order of magnitude. 



q 
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Fig. 14. — : The color gradient depicts the dust opacity as a function of the maximum grain size amax and 
the slope of the grain size distribution q. Results are shown for three different dust models (A, B and C) 
as described in the text. The minimum grain size is fixed at 5 x lO^'' cm. Thin solid curves show the dust 
opacity contours and are spaced by 0.5 dex. The thick solid curves show the possible pairs of {amax, q) that 
lead to values for the opacity slope /? equal to 0.5 (DG Tau) and 0.7 (RY Tan). The letters H and L indicate 
the high and low dust opacity models adopted in fitting the observed dust emission. 
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